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Abstract 

Using the Hugenholtz-Van Hove theorem, we derive general expressions for the 
quadratic and quartic symmetry energies in terms of single-nucleon potentials in 
isospin asymmetric nuclear matter. These analytical relations are useful for gaining 
deeper insights into the microscopic origins of the uncertainties in our knowledge 
on nuclear symmetry energies especially at supra-saturation densities. As examples, 
the formalism is applied to two model single-nucleon potentials that are widely used 
in transport model simulations of heavy-ion reactions. 
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1 Introduction 



One of the central issues currently under intense investigation in both nu- 
clear physics and astrophysics is the Equation of State (EOS) of neutron- 
rich nuclear matter [Tp|3] . For cold nuclear matter of isospin asymmetry 
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S = (Pn — Pp)/{Pn + Pp) at density p, the energy per nucleon E{p,S) can be 
expressed as an even series of 6 that respects the charge symmetry of strong 

interactions, namely, E{p,6) = Eo{p,0) + J2i=2Afi... EsymAp)^' where Esym,i{p) 
is the so-called symmetry energy of the ith order [3] and Eq{p, 0) is the EOS 
of symmetric nuclear matter. The quadratic term Esym,2{p) is most important 
and its value at normal nuclear matter density po is known to be around 
30 MeV from analyzing nuclear masses within liquid-drop models. Essen- 
tially, all microscopic many-body calculations have indicated that the higher- 
order terms are usually negligible around po, leading to the so-called empirical 
parabolic law of EOS even for 6 approaching unity for pure neutron matter. 
The Esym,2{p) is then generally regarded as the symmetry energy. For instance, 
the value of the quartic term has been estimated to be less than 1 MeV at po 
[H[5] . However, the presence of higher-order terms at supra- saturation densities 
can significantly modify the proton fraction in neutron stars at /3-equilibrium 
and thus the cooling mechanism of proto-neutron stars [5117] . It was also found 
that a tiny quartic term can cause a big change in the calculated core-crust 
transition density in neutron stars [39]. Therefore, precise evaluations of the 
quartic symmetry energy in neutron-rich matter are useful. Although much 
information about the EOS of symmetric nuclear matter Eq{p,0) has been 
accumulated over the past four decades, our knowledge about the density 
dependence of Esym,i{,p) is unfortunately still very poor. It has been gener- 
ally recognized that the Esym,i{,p), especially the quadratic and quartic terms, 
is critical for understanding not only the structure of rare isotopes and the 
reaction mechanism of heavy-ion collisions, but also many interesting issues 

mine the Esym,i{p) in neutron- rich matter has recently become a major goal 
in both nuclear physics and astrophysics. While significant progress has been 
made recently in constraining the Esym,2{p) especially around and below the 
saturation density, see, e.g., pnil8|19|l2D] . much more work needs to be done 
to constrain more tightly the Esym,i{p) at supra-saturation densities where 

dedicated experiments using advanced new detectors have now been planned 
to investigate the high density behavior of Esym,2{p) at many radioactive beam 
facilities around the world, it has become an urgent task to investigate theo- 
retically more deeply the fundamental origin of the extremely uncertain high 
density behavior of -Esym,2(p)- It is also of great interest to evaluate possible 
corrections due to the EsymA{p) term to the equation of state of asymmetric 
nuclear matter. 

Among existing proposals for extracting information about Esym,i{p) using ter- 
restrial laboratory experiments, transport model simulations have shown that 
many observables in heavy-ion reactions are particularly useful for studying 
Esym,i{p) in a broad density range. In these transport model simulations of 
heavy-ion reactions, the EOS enters the reaction dynamics and affects the 
final observables through the single- nucleon potential Un/p{p,S,k) where k 
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is the nucleon momentum. Except in situations where statistical equiUbrium 
is estabhshed and thus many observables are directly related to the binding 
energy E{p, 6) after correcting for finite-size effects, what is being directly 
probed in heavy-ion reactions is the single-nucleon potential f/„/p(p, 5, k). The 
latter is, however, directly related to the sjmimetry energy Esym,2{p) through 
the underlying nuclear effective interaction as first pointed out by Brueckner, 
Dabrowski and Haensel [35|36] using K-matrices within the Brueckner theory. 
They showed that if one expands Un/p{p, S, k) to the leading order in 6 as in 
the well-known Lane potential [57], i.e.. 



where Uo{p,k) and Usym,i{p,k) are, respectively, the nucleon isoscalar and 
isovector (symmetry) potentials, the quadratic symmetry energy is then [35)136] 




where t(/c_F) is the nucleon kinetic energy at the Fermi momentum kp = 
(37r^p/2)^/^ in symmetric nuclear matter of density p. The above equation 
indicates that the symmetry energy Esym,2{p) depends only on the single- 
particle kinetic and potential energies at the Fermi momentum kp- This is not 
surprising since the microscopic origin of the symmetry energy is the difference 
in the Fermi surfaces of neutrons and protons. The first term E'^^ = = 

is the trivial kinetic contribution due to the different Fermi mo- 
menta of neutrons and protons; the second term -kp is due to the 
momentum dependence of the isoscalar potential and also the fact that neu- 
trons and protons have different Fermi momenta; while the term \Usym{Pi kp) 
is due to the explicit isospin dependence of the nuclear strong interaction. 
For the isoscalar potential f/o(p, /c), reliable information about its density and 
momentum dependence has already been obtained from high energy heavy- 
ion collisions, see, e.g., ref. [13], albeit there are still some rooms for further 
improvements, particularly at high momenta/densities. On the contrary, the 
isovector potential Usym,i{p, k) is still not very well determined, especially at 
high densities and momenta, and has been identified as the key quantity re- 
sponsible for the uncertain high density behavior of the symmetry energy as 
stressed in ref. [3]. 

In the present work, we first show, using both the differential and integral 
formulations of the Hugenholtz-Van Hove (HVH) theorem [38], that the rela- 
tion in Eq.([5]) is valid in general. We then derive an expression for the quar- 
tic symmetry energy Esym,A{.p) in terms of the single-nucleon potential by 
keeping higher-order terms in the expansion of both the EOS and the single- 
nucleon potential. Applying the HVH formalism to two model single-nucleon 
potentials, namely, the Bombaci-Gale-Bertsch-Das Gupta (BGBD) potential 
[TB] and a modified Gogny Momentum-Dependent-Interaction (MDI) [251^0] . 



Un/p{p, 5, k) ^ Uo{p, k) ± UsyniAPi 



(1) 
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which are among the most widely used ones in studying isospin physics based 
on transport model simulations of heavy-ion reactions [SfT^ . we examine the 
relative contributions from the kinetic and various potential terms to Esym,2{p) 
and Esym,i{p)- We put the emphasis on identifying those terms that dominate 
the high density behaviors of Esym,2{p)- Finally, we evaluate the relative im- 
portance of the Esym,i{p) term by studying the Esym,i{,p) / Esym,2{,p) ratio as a 
function of density. 

The paper is organized as follows. In Section [2], based on the HVH theorem 
we derive general expressions for the higher-order symmetry energy terms 
Esym,2{p) and Esym,i{p) in terms of the single-nucleon isoscalar and isovector 
potentials. The derivation is carried out in Section 12.11 using the differential 
form of the HVH theorem by starting from the neutron and proton chemical 
potentials and in Section 12.21 using the integral form of the HVH theorem by 
starting from the total energy density of the system. Numerical results and 
discussions for both the BGBD and MDI interactions are given in Sections 13. II 
and 13.21 respectively. Finally, we give a summary in Section SJ 



2 Symmetry energy in terms of the single-nucleon potential 



In this section, we present two alternative approaches to derive expressions for 
the quadratic and quartic symmetry energy terms Esym,2{p) and Esym,A{,p)- 
Both are based on the Fermi gas model of interacting nucleons and satisfy 
the HVH theorem that was first derived in Ref. [5B]. There are, however, 
some technical differences between the two approaches in that a Taylor-series 
expansion is made on the single-nucleon energy in one approach but on the 
total energy of the system in the other. 



2.1 Derivation using the Hugenholtz-Van Hove theorem 



The Hugenholtz-Van Hove theorem [38] describes a fundamental relation among 
the Fermi energy Ep, the average energy per particle E and the pressure of 
the system P at the absolute temperature of zero. For a one-component sys- 
tem, in terms of the energy density = pE, the general HVH theorem can be 
written as [38)139] 



dipE) ^ dE ^ ^, 
Ef = ^ = =E + p— = E + P/p. (3) 

dp dp dp 

The above relation has been strictly proven to be valid for any interacting self- 
bound infinite Fermi system. It does not depend upon the precise nature of the 
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interaction. In the special case of nuclear matter at saturation density where 
the pressure P vanishes, the average energy per nucleon becomes equal to the 
Fermi energy, i.e., Ep = E. It is worthwhile to stress that the general HVH 
theorem of Eq.([n]) is valid at any arbitrary density as long as the temperature 
remains zero [38|39] . In fact, a successful theory for nuclear matter is required 
not only to describe satisfactorily all saturation properties of nuclear matter 
but also to fulfill the general HVH theorem at any density. In the following, 
we use the general HVH theorem to derive the relation between the nuclear 
symmetry energy and the single-nucleon potential. 

According to the HVH theorem, the chemical potentials of neutrons and pro- 
tons in isospin asymmetric nuclear matter of energy density C,{p, S) = pE{p, 6) 
are, respectively [38|39] . 



m) + u^ipj,k-F) = -^, (4) 

OPn 

t{e^) + u,{p,6,ep) = ^, (5) 

OPp 

where t{k) = hk'^/2m is the kinetic energy and Un/p is the neutron/proton 
single-particle potential. The Fermi momenta of neutrons and protons are 
kp = kpil + Sy/^ and kp = kp^l - Sy^'^, respectively. Subtracting Eq.® 
from Eq.(jl]) gives 



im) - m)] + [Unip, 5, k^p) - Up{p, 5, Fp)] = |l - |l. (6) 

(^Pn ^Pp 

The nucleon single-particle potentials can be expanded as a power series of 
5 while respecting the charge symmetry of nuclear interactions under the ex- 
change of neutrons and protons. 



Un{p,5,k) = Uo{p,k) + ^ Usym,i{p,k)5'' 

j=l,2,3... 

= f/o(p, k) + Usym,l{P^ k)6 + Usym,2{k)6^ + •.• (7) 



Up{p,5,k) = Uo{p,k)+ Usym,i{p,k){-6y 

i=l,2,3... 

= Uo{p, k) - Usym,l{p, k)6 + Usym,2{p, k)6'^ - .... (8) 

If one neglects the higher-order terms (5^, S^,...), Eq.Q and Eq.® reduce 
to the Lane potential in Eq.([T]). Expanding both the kinetic and potential 
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energies around the Fermi momentum kp, the left side of Eq.(|6]) can be further 
written as 



X [( E Fuwr-i E mi-syy] 

i=l,2,3.. i=l,2,3.. 
+ E U,ymAp^kF)[6'-i-6y] 

1=1,2,3... 

, 1 d'Usym,l{p, w 

/=1,2,3... j=l,2,3... ^'^ 

X [( E E mi-syn-sy] 

jr=l,2,3.. j=l,2,3.. 

.2a[t(fc) + f/o(p, k)] 



'3 dk 



kF + 2Usy„^,lip,kF)]6 + ..., (9) 



where we have introduced the function F{j) = jt[|(| — — j + !)]• For 

the right side of Eq.([6]), expanding in powers of 6 gives 



dpn dpp pd6 

= ^Esym,2{p)5 + 8E,j,„,4(p)5^ + 12E,y„,6(p)5^ + .... (10) 



Comparing the coefficient of each 5* term in Eq.(l9]) with that in Eq.f llOp then 
gives the symmetry energy of any order. For instance, the quadratic term 



, , ld[t{k) + Uo{p,k)]. , 1 ^ , , , 

Esym,2[P) = IkpKp + -Usym,l[P, Kp) 

= ^tikp) + Ifc^ -kp + ^Usym,l{p, kp) (11) 



is identical to that in Eq.(l2]), while the quartic term can be written as 



6 



E. 



sym,i 



5 d[t{k) + Uo{p,, 



324 



dk 



1 d'[t{k) + Uo{p,k)] 



108 



dk^ 



2 1 d^m + Uo{p,k)] ^ 3 



648 



1 dUsym,l{p,k) 



36 



1 d^U. 



sym,l 



72 



1 dUsymAP^k) 



12 



9A; 



sym,3 



dk^ 
(p, 



(12) 



^.^ Derivation using the total energy of an interacting Fermi gas 



The symmetry energy of any order obtained in the previous subsection can 
also be derived within the interacting Fermi gas model [^TM2] . Although the 
derivation is more tedious, it is physically interesting and mathematically in- 
structional. 

There are different ways to calculate the total energy of a Fermi system at a 
given density p. As explained in detail by Bertsch and Das Gupta |13], one 
could determine the total energy starting with empty space, adding particles 
until the desired density is reached. Each added particle would contribute an 
energy k{pxY /2m + U{px, k{px)), where k is the Fermi momentum correspond- 
ing to the density px of particles already added to the system. As stressed by 
Bertsch and Das Gupta, the single-particle potential U{px, k{px)) is not the 
same as the potential energy per particle as one might at first guess. In this 
way, the total energy density of the asymmetric nuclear matter written in 
coordinate space is 

^= / [k{Px) /2m + Un{px,k{px))]dpx+ / [k{px) /2m + Up{px, k{px))]dpx, 
Jo Jo 

(13) 

where p^ + fP = p. The single particle potential U{px, k{px)) for the neutron 
or proton in the first or second integral in Eq. (fT3l) can be rewritten as 



U{px,k{px)) = U{p,6*,k{px)), (14) 

where px = p(l + S*)/2 in terms of the local isospin asymmetry S* (to be 
discussed in detail in the following). The single particle potential can be further 
expanded as a series of 6* 



uip,6*,k{px)) = Uoip,k{px))+ J2 Usy^MKpMn (15) 

j=l,2,3... 

= [/o(p, k{px)) + UsymAp, k{px)){,S*) + Usyyn,2{p, k{px)){S*f + .... 
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Similar to the interpretation of Eq. (ll3p . one can raise the momenta of all 
particles in the system from zero to the Fermi momentum corresponding to 
the density p. In either the coordinate or momentum space, as one increases 
the density or momentum for both neutrons and protons to build up the 
desired Fermi system, the local isospin asymmetry 5* changes continuously as 
a function of density or momentum, i.e.. 



2p^ 
P 



^6* 



1,3 

A. 77 



(16) 



where is either the neutron or proton density used in the two terms of 
Eg. (1131) . and k is the neutron or proton momentum used in the two terms of 
Eq. ffTSl) below. At the respective Fermi surfaces of neutrons and protons, the 
local isospin asymmetry 6* reduce to the global one 









. p 












. p 







-5). 



(17) 



In order to make use of the same technique in expanding both the kinetic and 
potential terms in 6 series, it is more convenient to express the energy density 
in momentum space 



[t{k) + Un{p,S* 



Wdk + 



[t{k) + Up{p,6*,k)]k^dk 



(18) 



Before we proceed to derive expressions for Esym,2{p) and Esym,4{p), it is criti- 
cal to examine whether the energy density given in Eq. (fT3l) or Eq. ( |T8i) satisfies 
the HVH theorem. Noticing that 



dkp dp^' 



dk^^p 



I 



dk^'p^ 



dkl 



I 



37r2 



dk-^ 



dk'^p 



TX 



,(19) 



then it is straightforward to show using Eq. OlSp that 



ap" 



1 



-Am) + Ur,{p,6,kl)]{k 



n\2 



TT 



= t{kl) + Un{p.5, kl) 

Similarly, one can show that 



(20) 



8 



(21) 



Thus, the HVH theorem is indeed satisfied by the energy density expressed in 
both Eq.([T3]) and Eg.ffTS]). 



To obtain expressions for Esym,2{p) and Esym,i{p)i we make a Taylor-series 
expansion of the energy density in Eq. lITH]) . To proceed, it is useful to first 
recall that for any continuous function Y{k) the Taylor-series expansion of 

the integral f{k'^F) = /o^ Y{k)k dk around kp leads to [S] 



fik'p) = fikp) + Y{kp) ■ kl ■ {k^ - kp) 



dY 



kp-kl + 2kp ■ Y{kp) 



{k^-kpf + ..., 



(22) 



Moreover, it is easy to show that 



kp — kp 
kp — kp 



3 9 81 243 



3 9 81 243 



5' + ... 



(23) 



The kinetic energy per nucleon can then be expanded as 



T = a I ' t{k)k^dk + a / t{k)k^dk 
/o Jo 

rkp 



= 2a [ ^ t{k)k^dk + at{kp)kl\{k'l - kp) + {kl - kp)] 
Jo 

+ ^^^^^\kpm - kp)' + (Fp - kp)'] 



2! dk 

ad'[t{k)e] 

3! 



kplik^P - kp)' + (k^P - kp) 



+ ^^^^^^\kpm - kp)' + {Fp - kp)'] 



4! (9fc3 



1 dtik) 



2a I " t{k)edk+-^^ U„ kp6' 



5 dt{k) 
324 dk 



kp Kp 



kp — 



6 dk 
1 dHik) 



108 9P 



\kp i^p 



1 dH{k) 
648 dk^ 



\kp i^p 



s' + 



(24) 



where a = 3/{2kp). Thus, the kinetic contribution to the Esym,2{,p) is 
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E. 



(p) = -— rr^lt./c 



sym,2 



6 



-t{ki 



2\ 2/3 



.2/3 



6m \ 2 



(25) 



and its contribution to Esym,4{p) is 



rpkin / \ ^ dt{Ji) , 
1 



1 



108 9fc2 



81 ^ ~ 162m \ 2 j ' 



\kp i^F 



„2/3 



1 dH{k) ^ 
648 9A;3 



(26) 



For the potential energy per nucleon written in momentum space, we first 
expand the single-nucleon potential Un/p in local isospin asymmetry 5*{k) 



U = a Un{p,S\k)k^dk + a / UJp,6\k)k^dk 
Jo Jo 



rk^ 

Uo{p,k)k'^dk + a / Uo{p,k)k^dk] 
Jo 



+ E 



=1,2,3., 



a UsymAP^^WY^dk + a Usym4p,k){6*yk'dk 
Jo Jo 



k^ 



(27) 



We then expand each integral in the above equation around the Fermi mo- 
mentum kp using Eq. (l22l) . The Uq contribution to U is 



Ef;°\p) = a Uo{p, k)k^dk + a [''^ f/o(p, k)k'^dk 
Jo Jo 



2a rU^MkHk+^-^'^'^P^''^ 
Jo 



+ 



5 dU^{p,k) 
324 dk 



kp 



6 dk 
1 d'Uoip,k) 
108 dk^ 



kp5'^ 



\, k^ 



1 d'Uojp.k) 
648 dk^ 



, k^ 



6' + 



{2i 



Contributions from the higher-order potential terms in 5 can be obtained 
similarly. For the second order-symmetry energy Esym,2{p), only the Usym,i 
term contributes 
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k 



*'-^dk 
k 



■a / U,y„^Ap,k)[{—f-l]k^dk + a / Usym,i{p,k)[{—f-l]k^dk 

Jo Kp Jo Kp 

■-2a ['^ UsymAp,k)[{^f-l]edk 

Jo ftp 



+ -Usym,l{p, kF)S'^ + 



1 dUsym,i{p,k) ,1 d'^Usym,i{p,k) I 2 



36 dk 



72 



(29) 



To obtain the fourth-order symmetry energy Esym,i{p) , ah Usym,i, Usym,2, and 
Usym,3 terms are needed, i.e. 



ET{p) = aj^ UsymM kWYk'dk + aj^ f/s,^,2(p, k){6*yk'dk 

2 



a / Usym,2{p,k) 



2a / Usyni,2{p,k) 

•J 



kp) 



k 



k'dk + a / Usyn^AP^ k) [{—Y - Ifk'dk 

Jo ftp 



12 ok 



(30) 



and 



El°\p) = aj^ UsymM^ kWfk'dk + aj^ UsymAp^ k){d*yk'dk 

3 



a / Usym,?,{p,k) 



2a / Usym,z{p,k) 

J 



k \ 



' k \ 



hp 



k 



k'dk + a / UsymAP^ [{—Y - lYk'dk 

Jo Kp 



k'^dk + -t/sym,3(p, kp)5'^ + .... 



(31) 



Combining all coefficients of the 5^ and 5^ terms in both the kinetic and 
potential parts given above, one sees that the Esym,2{p) and Esym,4:{p) obtained 
here are exactly the same as those given in using directly the HVH theorem. 
Moreover, they agree with earlier results obtained by Brueckner et al. [35|l36] . 
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3 Applications and discussions 



As shown in the previous section, the symmetry energy can be explicitly 
separated into the kinetic energy term T and the potential terms Uq and 
Usym,i at the Fermi momentum kp. To evaluate their relative contributions 
to the symmetry energies, especially for the second-order and fourth-order 
terms Esym,2{p) and EsymA{p), we consider in this section two typical single- 
nucleon potentials that have been widely used in tansport model simulations 
of heavy-ion reactions. 



3.1 The Bombaci-Gale-Bertsch-Das Gupta potential 



As a first example, we use the phenomenological potential of Bombaci-Gale- 
Bertsch-Das Gupta [16] 



Ur{u, 6, k)=Au + Bu" - -(a - + xs)^ 

3 (J + 1 2 



± 



2 .,1 ^ 4 5 1 
'3A(- + xo)n---^(- + X3K 



5po L2 



5 



C±'^-^5]u.g{k), (32) 



where u = p/po is the reduced density and ± is for neutrons/protons. In 
the above, we have I, = [2/{27rf] J d^kfr{k)g{k) with g{k) = 1/[1 + (|)2] 
being a momentum regulator and /r(/c) being the phase space distribution 
function. The parameter A has the value A = 1.5kp, where kp is the nucleon 
Fermi wave number in symmetric nuclear matter at pQ. With A=-144 MeV, 
B=203.3 MeV, C=-75 MeV and a = 7/6, the BGBD potential reproduces 
all ground state properties including an incompressibility -ft'o=210 MeV for 
symmetric nuclear matter [16]. The three parameters Xo,X3 and Zi can be 
adjusted to give different symmetry energy Esym,2{p) and the neutron-proton 
effective mass splitting m* — m* [T6]l44|45j . For example, the parameter set 
zi = —36.75 MeV, xo = —1.477 and x^ = —1.01 leads to m* > m* while the 
one with zi = 50 MeV, xo = 1.589 and x^ = —0.195 leads to m* < m* at all 
non-zero densities and isospin asymmetries. 



On expanding the BGBD potential in S, the coefficients of the first few terms 
are 
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[/o(p, k) = Un/p\ 5=0 

= Au + Bu" 



2CA2 

PO 



kp — Atan 



' kp^ 

T 



Usym,l{P^ k) = ± 



1 9f/n/p I 



1! 86 

-W- 

3 ^2 

C-8zi 



5=0 



Xo)u 

u ■ gi^kp) + 



4 B 

3a 

C- 



- r2 

8^1 



u-g{k), 



Usym,2{p, k) 



1 d^U, 



n/p I 



2! 



1) 



5=0 
B 



-A 



J-^9ikp) 



Usym,3iP, k) 



1 d^U, 



n/p I 



3! 
C-8 

135 



<5=0 



1,2 rri,2 
i^p / ^Ji^p 



Thus, the second-order symmetry energy Esym,2{p) is given by 



Esym,2{p) = ^t{kp) + IfcF -f^F + ^Usyra,l{p, kp) 



/37r 
6m ' 



2\ 2/3 



P 



2/3 



1 ./I 



- -^^^9{kF) 

2 B A 
+ X3 n"" 



C-8zi 



u ■ gikp), 



and the fourth-order symmetry energy EsymA^p) is 
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0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



P/Po 

Fig. 1. The kinetic energy part (T), the isoscalar potential part (C/q) and the isovec- 
tor potential part {Usym,i) of the symmetry energy Esym,2 from the BGBD potential 
with m* > m* (left) and for m* < m* (right). 

In Fig. [T]we compare Esym,2{p) and its three components in the two cases of 
m* > m* and m* < m*. It is seen that the kinetic and isoscalar contribu- 
tions are the same in both cases. However, they have significantly different 
isovector potentials Usym,i, leading thus to different Esym,2{,p) especially at 
supra-saturation densities. 

Various contributions to the fourth-order symmetry energies Esym,i{p) in the 
two cases are compared in Fig. [21 Similar to Esym,2{p), the contributions of 
the T and Uq terms to Esym^4^{p) are positive and they are the same in both 
cases. Interestingly, the Usym,i term also plays the most important role in 
determining the high-density behavior of Esym,ii.P)- It is positive in the case 
of m* > m* but negative in the case of m* < m*, resulting in very different 
behaviors oi Esym,A{p) at supra-saturation densities. Moreover, it is interesting 
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Fig. 2. The kinetic and various potential contributions to the fourth-order symmetry 
energy Egym^A, with the BGBD potential for ?ti* > m* (left) and for < ?ti* (right). 



to note that Esym,4:{p) receives no contribution from the Usym,2 term. This is 
not surprising because the Usym,2 term in the BGBD interaction is momentum 
independent and its contribution to Esym,i{p) is actually ■f^ '^^^'^aif^''^ ^^=0- 
On the contrary, the Usym,3 term still contributes to EsymA{p) via \Usym,3{kp) 
although it is momentum independent too. In the two cases considered here, 
the contributions from the Usym,3 term also have opposite sign. 

To compare the fourth-order term Esym,A{p) with the second-order term Esym,2{p) 
more clearly, we show in Fig. [3] their ratio Esym,4:{p) / Esym,2{p) as a function 
of the reduced density p/po- Obviously, the relative value of Esym,4:{p) is gen- 
erally small. However, it can reach up to about ±10% at high densities for 
both cases of m* > m* and m* < m* . It may thus lead to an appreciable 
modification in the proton fraction and therefore the properties of neutron 
stars at /^-equilibrium. 



3.2 A modified Gogny Momentum-Dependent-Interaction 



In this subsection, we discuss the symmetry energy obtained from the MDI 
interaction [23], which is derived from the Hartree-Fock approximation using 
a modified Gongy effective interaction |1D] 
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Fig. 3. The ratio of Esym,A over Esym,2 with the BGBD potential for m* > m* (left) 
and for m* < m* (right). 



U{p,6,p,T) = Au{x)^ + A{x)^ 
Po Po 

Po cr 



T^^Spr' 

1 Po 



2a 



Po 



1 + (P-^)VA^ 



2Cr, 

Po 



d^p' 



fr'{r,f) 



1 + [p-p'Y/A^' 



(39) 



In the above, r = 1/2 (—1/2) for neutrons (protons) and tt^t' ;(T = 4/3 
is the density- dependence parameter; frir^p) is the phase space distribution 
function at coordinate r and momentum p. The parameters Cr,T, Cr,T' and 
A are obtained by fitting the nuclear matter saturation properties [23]. The 
momentum dependence of the symmetry potential stems from the different 
interaction strength parameters Cry and Cr^r for a nucleon of isospin r in- 
teracting, respectively, with unlike and like nucleons in the background fields. 
More specifically, CunUke = —103.4 MeV while Cuke = —11.7 MeV. The quan- 
tities Au{x) = —95.98 — x-^^ and Ai(x) = —120.57 + are parameters. 
The parameters B and a in the MDI single-particle potential are related to 
the to and a in the Gogny effective interaction via to = I ^ and a = a + 1 
[iO] . The parameter x is related to the spin(isospin)-dependence parameter xq 
via X = {1 + 2xo)/3 [16]. On expanding the single-nucleon potential in 6, the 
first four terms are 
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According to Eq. lfTTl) . the second-order symmetry energy ii/sym,2(p) is 



17 



■kp + -Usym,i{p, kp) 



2\ 2/3 



P 



2/3 



+ 



6m Y 2 
3po 



3po 



/l3 



A2 



(44) 



and according to Eq.f fT2|) the fourth-order symmetry energy -£'gy^ 4(p) is 
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(45) 



As one expects, the above expressions are identical to those derived directly 
from the exact MDI EOS using [17] 



£'sym,2lPj - ^ ^ 1 5=0 

^sym,4(Pj - ^ |5=0- 



(46) 



In Fig. m we show the kinetic (T), isoscalar (Uo) and isovector {Usym,i) poten- 
tial contributions to Esym,2 for the three different spin (isospin)-dependence 
parameter a; = 1, 0, and -1. We notice that the kinetic (T) and the isoscalar 
potential (Uo) contributions are the same for the three different x values. 
As pointed out in Ref. [23], it is the isovector potential Usym,i that is caus- 
ing the different density dependence of Esym,2- For instance, with x = 1 the 
Usym,i term decreases very quickly with increasing density and thus results in 
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Fig. 4. The kinetic energy part (T), the isoscalar potential part {Uq) and the isovec- 
tor potential part {Usym,i) of the symmetry energy Esym,2 from the MDI interaction 
with X = 1, and -1. 

a super-soft symmetry energy at supra-saturation densities. On the contrary, 
the symmetry energy Esym,2 at supra-saturation densities is very stiff for both 
X = and x = — 1 as the contribution of the Usym,i term becomes very positive 
with smaller values of x. 
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Fig. 5. The kinetic energy and potential contributions to the fourth-order symmetry 
energy Esym,4 from the MDI interaction. 
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Unlike the second-order term Esym,2, the fourth-order symmetry energy Esym,^ 
is independent of the spin (isospin)-dependence parameter x. Shown in Fig. 
E] are the various contributions to the fourth-order symmetry energy Egym^i- 
Comparing these with the results obtained using the BGBD in Fig. [21 we find 
that the T and Uq terms from these two interactions are almost identical. 
However, there exists some differences for other terms. For the MDI interac- 
tion, the Usym,2 term is negative and becomes very important for determining 
Esym,i- One the contrary, the contributions from the Usym,i and Usym,3 terms 
are positive and they are relatively small as compared to Usym,2- Generally, 
the behavior of Esy„^^4^ from the MDI interaction is very similar to that from 
the BGBD interaction for the case of m* > m*. 
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Fig. 6. The ratio of Esym,4: over Esym,2 obtained from the MDI interaction as a 
function of reduced density p/po for x = 1, 0, and -1. 

To compare more directly Esym,i with Esym,2, their ratio Esym,^/ Esym,2 is plot- 
ted in Fig. [6] as a function of reduced density for x = 1, 0, and -1. It is seen 
that with X = 1 there is a sharp break in the curve around 3po- This is because 
the second-order symmetry energy Esym,2 changes from positive to negative 
around 3po in this case. However, this is not the case for both x = —1 and 
a; = where Esym,2 remains positive at all densities. In all cases, -Esym,4 is very 
small compared to Esym,2- For both BGBD and MDI interactions, the small 
values of Egym^i up to several times the normal density clearly shows that 
the parabolic approximation of the EOS is well justified for most purposes. 
However, cares have to be taken in evaluating the core-crust transition density 
where the energy curvatures are involved [9]. 
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4 Summary 



In summary, using the Hugenholtz-Van Hove theorem we have derived gen- 
eral expressions for the quadratic and quartic symmetry energies in terms 
of single-particle potentials in isospin asymmetric nuclear matter. Identical 
results are obtained by using two approaches, i.e., one based on the single- 
particle potential and the other based on the total energy, that are physically 
identical although mathematically different. By using the derived analytical 
formulas, the symmetry energies are explicitly separated into the kinetic and 
several potential parts. The formalism is applied to two typical single-nucleon 
potentials, namely the Bombaci-Gale-Bertsch-Das Gupta (BGBD) potential 
and the modified Gogny Momentum-Dependent-Interaction (MDI), that are 
widely used in transport model simulations of heavy-ion reactions. We find 
that for both interactions the isovector potential is responsible for the un- 
certain high density behavior of the quadratic symmetry energy. Also, the 
magnitude of the quartic symmetry energy in both cases is found to be sig- 
nificantly smaller than that of the quadratic symmetry energy. We expect 
that the analytical formulas for the nuclear symmetry energies derived in the 
present study will be useful in extracting reliable information about the EOS 
of neutron-rich nuclear matter from heavy-ion reactions. 
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